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ABSTRACT 

Moments of the BBGKY equations for spatial correlation functions of 
cosmological density perturbations are used to obtain a differential equation 
for the evolution of the dimensionless function, h = —(v/ax), where v 
is the mean relative pair velocity. The BBGKY equations are closed 
using a hierarchical scaling ansatz for the 3-point correlation function. 
Scale-invariant solutions derived earlier by Davis and Peebles are then used in 
the non-linear regime, along with the generalised stable clustering hypothesis 
[h — > const.), to obtain an expression for the asymptotic value of h, in terms 
of the power law index of clustering, 7, and the tangential and radial velocity 
dispersions. The Davis-Peebles solution is found to require that tangential 
dispersions are larger than radial ones, in the non-linear regime; this can be 
understood on physical grounds. Finally, stability analysis of the solution 
demonstrates that the allowed asymptotic values of h, consistent with the 
stable clustering hypothesis, lie in the range < h < 1/2. Thus, if the 
Davis-Peebles scale-invariant solution (and the hierarchical model for the 
3-pt function) is correct, the standard stable clustering picture (h — > 1 as 
£ — > 00) is not allowed in the non-linear regime of structure formation. 



Subject headings: cosmology: theory — large scale structure of the Universe 



1. Introduction 



The cosmological BBGKY equations have been used on a number of occasions to 
study the evolution of non-linear density fluctuations in an expanding, flat (Q = 1), 



[background Universe (Pavis fc Peebles 1977 : 



Huamsuwan fc Fry 1992|; 



Yano fc Goudaj 
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1997). These equations have the advantage of dealing directly with statistical quantities, 
i.e. the spatial N-point correlation functions; however, they form an infinite hierarchy 
in which the equation for the N-point correlation function contains terms involving the 
(N+l)-point function. This necessitates the development of closure schemes in which 
the hierarchy is cut off at a finite number of equations by making some assumptions 
wherein the higher order correlation functions are written in terms of the lower order 
ones. Such a scheme was used by Davis & Peebles (1977; hereafter DP) to demonstrate 
the existence of a similarity solution to the equations; the assumptions used were the 
stable clustering hypothesis, the vanishing of velocity skewness, a hierarchical model for 
the 3-point correlation function, ( oc< £ 2 >, where £ is the 2-point correlation function, 
and finally, a factorisation of the 2- and 3-body phase space distributions which gives a 
form for the 3-body weighted pair velocity. In this solution, an initial power spectrum 
of form P Q {k) ~ k n evolves to give a 2-point correlation function 



eWocr- 7 = ^7^ (1) 
(n + 5) 

in the strong clustering regime with £ 3> 1. Ruamsuwan & Fry (1992; hereafter RF) 
showed that the assumption of vanishing velocity skewness was not a priori necessary 
but could be derived as a result of other, more general assumptions; the solution was 
also shown to be marginally stable to perturbations. Finally, Yano & Gouda (1997) 
found that the stability condition used by DP was satisfied for the unique case of 
vanishing skewness and the specific form for ( in terms of the products of the 2-point 
correlation functions. The power index of the 2-point function in the strong clustering 
regime depends, in general, on the mean relative physical velocity, the skewness and the 
3-point correlation function. 

The present work considers the evolution of the function, h = — (< v >/ax), i.e. 
the ratio of the mean relative peculiar velocity to the Hubble velocity. The BBGKY 
equations are not derived again; instead, we use the relevant moment equations from RF 
to proceed. The form for the 3-point correlation function is the same as in DP and RF. 
We derive an equation for the evolution of h in terms of the mean 2-point correlation 
function, £, using the ansatz that h is a function of £ alone ( |Hamilton et al. 1991| ; 



Nityananda fc Padmanabhan 1994| , hereafter NP). This ansatz is, however, not crucial 



to the later discussion. (Note, further, that the above equation is derived using only the 
zeroth and first moments of the second BBGKY equation and hence does not contain 
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any assumptions regarding the form of the velocity skewness or the higher velocity 
moments, excepting the fact that they be such as to yield the DP similarity solution 
for the two-point correlation function and the velocity dispersions; the assumption of 
vanishing velocity skewness is not used here. Of course, as shown by RF, this assumption 
is not necessary to arrive at the DP solution; in fact, zero skewness results as only one 
particular case in their closure scheme.) The DP solution is then substituted in the 
equation for h and the generalised stable clustering hypothesis (h — > const., for £ 3> 1) 
used, to obtain an expression relating the asymptotic value of h to 7, the power law 
index of clustering, and the radial and tangential pair velocity dispersions. This, and 
the requirement that h is real, gives the constraint that tangential dispersions are larger 
than radial dispersions, in the non-linear regime. Finally, a stability analysis carried out 
by perturbing about the above solution shows that stable solutions are attained only for 
asymptotic values of h in the range < h < 0.5. Thus, the standard stable clustering 
picture (h — > 1 as £ — > 00) is incompatible with the DP scaling solution. 



2. The BBGKY hierarchy 

The zeroth and first moments of the 2 nd BBGKY equation are respectively (RF, 
equations (22) and (23)) 



12 



1 d 



+ — 



dt a dx 



<v\ 2 >(l + C. 



12 



12 



(2) 



ld_ 

a dt 



a < v\ 2 > + f 



12, 



+ - 



1 d 



a dx 1 



12 



< ^12^12 > (1 + 62) 



+ Gp b a I d 3 x 3 (£i 3 + 63 + Cm) ' ' ' ' ' 



Xl3 | 3 I X 2 3 | 3 



(3) 



where x 



12 



.1 I 



Xc 



denotes the separation between particles labelled 1 and 2 and 



v\ 2 = v\ — v\ is their relative velocity (i is a vector index). £12 and £123 are the usual 2- 
and 3-point correlation functions respectively while pb is the background density. 
The number of independent vector components in equation (|3p can be reduced by 
making use of various cosmological symmetries. The isotropy of the background 



- 4 - 



Universe implies that the mean relative pair velocity, < v\ 2 >, should be aligned with 
the pair separation, i. e. 



< v\ 2 >= v(x 12 )x\ 2 



(4) 



where x % denotes a unit vector along x l , and the pair velocity dispersion should have 
longitudinal and transverse polarisations about the mean II and S (RF). 



< Av i Av j >= + (5 ij - x l x 3 )E(x) 



(5) 



We note that II and S are peculiar velocity dispersions. Finally, C123 is taken to be the 
hierarchical form (DP; RF) 



Cl23 — Q (£l2£l3 + ^13^23 + ^126 



23/ 



(6) 



This form for (123 satisfies the necessary symmetry under the exchange of indices; 
further, it vanishes when any of the three points is removed to a large distance. Using 
the above ansatz, we obtain 



d£ 1 d r 
— H 

dt ax 2 dx 



x 2 v(l + 



(7) 



and 



ld_ 

a dt 



au(l + 



+ 



1 r x 2 (n + t; 2 )(i + ol--s(i + o 

J ax 



+ 



ax 2 dx 
2Gpba 



x* Jo 



x p COS 

dh£{z) + 2GQ Pb a / d 3 z £(x) + £(z) f (z - x)— 5- = (8) 

J z 



We next differentiate equation ([?]) with respect to t and equation (|8]) with respect to x; 
the resulting equations can be combined to yield 



d r 2 d(x 2 ^)i d 



dt 



dt 



+ 



dx 2 



x 2 (n + ^ 2 )(i + o 



<9x 



2Ex(l + 



- 5 - 



+ 



d_ 

dx 



2G Pb a 2 f ' d?z£(z) 
Jo 



+ 



d_ 

dx 



2GQp b a 2 x 2 M =0 (9) 



where M is defined by 



M — d 3 z 



+ / \ xi x cos 



(10) 



Now, the mean 2-point correlation function, a), is defined by 



- 3 r x 

£(x, a) = — / dx£(x, a)x 2 



x 3 Jo 



(11) 



Substituting for £ in terms of £, we obtain 



d_ 
~dt 



2 d ( d fX 3 C\\ 

a dt\d x y 3 n 



+ 



dx 2 



<9x 



x 3 (l + 



9a; 



2E 9 rx 3 (l +0 
a; <9a? 



+ 



d_ 

dx 



8nG 



d 



p b a 2 x% + 2GQp b a 2 x 2 M = (12) 
J ox 1 " J 



After integrating over x and carrying out some algebra, this gives 



dA 2 Id A 



3?-(fc| + fc 2 ) 



<9F 



<9X 



9 



4F+— — F——[hj + h \ + 2h ± F — 



dX 



9QMe 



-x 



(13) 



In the above, we have defined F by 



F = *g + 3(l+f) 



(14) 



, substituted X = In a; and A = In a, and defined 
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However, d£/dA = hF (NP). Hence, 



dh 



dh 



hF 



OF- 



dhl 



dx- 



dX 



9QMe~ x 
Aix 



(16) 



Here, we make the ansatz, h = ( Hamilton et al. 1991| ; NP; |Mo, Jain fc White 1995 
Padmanabhan fc Engineer 1998| ). This gives 



dh 



dA dX 



h , 2 



«91nF 
<9X 



9MQe 
AttF 



-x 



(17) 



and, finally, 



h l2 3£ , 2/ <91nF x dhf, _ 9MQe" 
where we have used (NP) 



H_ ft ^ = 3 M l + ?) (19) 

Further, the ansatz, h = has been used to convert the partial derivative into a 
total derivative. Equation ( |I8D governs the evolution of ft, in terms of the mean 2-point 
correlation function £. 



3. The non-linear regime 

We will consider the small separation, strong clustering limit, £ 3> 1. In this 
regime, the 2-point correlation function was found by DP to exhibit a scale- invariant, 
power-law behaviour, with the assumption of stable clustering. The latter is physically 
well-motivated as it seems reasonable to expect stable, bound systems to form under the 
influence of gravity. Such systems would neither expand nor contract and would hence 
have peculiar velocities equal and opposite to the Hubble expansion, i.e. v l pec = —ax 1 . 
The stable clustering ansatz has, however, not been deduced from any fundamental 
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considerations and might certainly be considered suspect if mergers of structures are 
important ( |Padmanabhan et al. 1996| ). Also, while N-body simulations ( [Hamilton et 
al. 1991 ) indicate that h — > 1 for £ ^> 1, the results are, at best, inconclusive. The 
argument can, however, be generalised ( Padmanabhan 1997 ) as the function h, which is 
the ratio of two velocities, should tend to some constant value, if the virialised systems 
have reached stationarity in the statistical sense. We will use this more general stability 
condition, namely h — > const., to proceed; using this assumption, equation (0) reduces 
to 



<9£ h d 

da x 2 dx 



x s ti 



(20) 



The above equation has the power-law solution, £ ~ a^x 7 with f3 = (3 — j)h. 
Dimensional analysis of equation @ reveals that II and £ scale as II 
Now, for a flat (Q = 1) Universe, a 2 oc a -1 . This implies that 



n 

a 2 x 2 



oc a^-^x^ 



IU 



(21) 



and, similarly, 



hi 



(22) 



where n o and S D are constants of proportionality. Also, ( |Yano fc Gouda 19971 ) we can 
write M = M'xtf , where M' is another constant. In the non-linear limit, £ ^> 1, 



^=(3-7)^ + 3 



(23) 



and 



<91nF 

dX 



-7 



1 + 



5 



3-7 /v r 



(24) 
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Substituting for F, M, h\ and £ in equation (ITS]) and using the limit £ ^ 1, we obtain 



- /> 2 + | - (^-)D( 7 ) = ?C( 7 ) + (25) 
d£ 2 3-7 £ 



where I? (7) and C7( 7 ) are defined by 



and 



C( T ) = 2H (2 - 7) - 2S + Jg M '^ (27) 

Equation (|25|) is exact (within the exact power law solutions for £, /z|, and M) upto 
order constant, in the limit of large £, with terms of order $(l/£) neglected. 
Clearly, C( 7 ) must be exactly zero, as otherwise, h oc y£ for £ ^> 1, which violates the 
stable clustering hypothesis. Note that C( 7 ) ~ is not sufficient as, if this were the 
case, the term in C{j) would cause h to grow with £, for sufficiently large £. Thus 



2S °- 2n °( 2 -^ = 4^) (28) 
This equation is equivalent to equation (48) of RF, albeit in slightly different form. 



Equation (p5|) thus reduces to 



3^§ - h> + \ = A (29) 
d£ 2 

where A = [3/(3 - 7)]D( 7 ). 

Now, the ansatz of stable clustering implies that ^dh/dC, — > as £ — > 00. Equation (|29| ) 
then gives 



, 9 /i 

h — - + A = 



(30) 



- 9 - 



i.e. 



4 



1 ± Vl - 16A (31) 
Since ft, is a real quantity, the above equation immediately yields 



(32) 



16 

This gives 



9QM' /3- 7 \ 
1+7n °-47(3^)^l^) (33) 



Replacing for 9QM'/47r(3 — 7) from equation (|28|), we obtain 



2£ >(4- 7 )n + l-( 3 - I ^) (34) 

RF have pointed out that various integrals in the BBGKY hierarchy do not converge 
unless < 7 < 2. This range of 7 values implies that 1 — (3 — 7)/48 > and (4 — 7) > 2. 
Thus, the inequality (53) gives 



s > n (35) 

Thus, the DP solution requires that tangential dispersions exceed radial dispersions 
in the non-linear regime. This is understandable on physical grounds, as tangential 
dispersions cause deviations from radial infall; stable structures would hence only be 
expected to form once these dispersions become comparable to or larger than the radial 
ones. 

The general solution to equation (^9|) has the following three different forms depending 
on whether (1 — 16 A) is positive, negative or zero. The constant of integration is denoted 
by B in each case. 



- 10 - 



1. For (1 - 16 A) < 0, the solution is 



In ( Bf l& 



In 



(2h 2 -h + 2A 



1/4 




2 V 7 + 45 

2. Next, for A = 1/16, the solution has the form 



/ 3-7 
7 + 45 



(36) 



2/3 



2/i 2 - /i + 2A 



exp 



3. Finally, for (1 — 16A) > 0, the solution is 



4/i-l 



(37) 



2h 2 - h + 2A 



p-Ah + 1 
p + 4h-l 



i/p 



(38) 



where we have defined p = y/T— 1QA. As mentioned earlier, (1 — 16 A) cannot be 
negative since h is a real quantity; equation can hence be immediately ruled out 
as a possible solution. Further, although £dh/d!; — > 0, h = 0.25(1 ± \p\) is certainly a 
solution of equation (|29|) , it is not clear if this can be actually reached from the general 
solution embodied in equations fl3"7| ) and (R3). We hence perturb these solutions by 
writing h = h + e, where e is the perturbation parameter and h Q satisfies the equation 



hi- — 
° 2 



7 



(39) 



i.e. 



h = 0.25 (1 ± \p\) 



(40) 



We then attempt to impose the condition e — * as £ — > oo; if this is possible, it 
clearly indicates that h — > /i as ^ — > oo, i.e. a solution exists which satisfies the stable 
clustering hypothesis. We initially consider the case A = 1/16 (p = 0) and rewrite 
equation (|3T|) as 



5£ 



■2/3 



2/t 2 - h + 2A + Ah e + e 2 



cxp 



1-Ah- 4e. 



(41) 
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Since h a = 0.25 for p = 0, this gives 



2/3 



e 2 " 




■ 1 ■ 




exp 


~2e. 



(42) 



One can impose the condition e — > as £ — > oo, in the above equation, only if e is 
negative. A = 1/16 is thus an allowed solution and h — > 0.25 from below, as £ — > oo. 
Finally, we consider the case (1 — 16 A) > 0, i.e. A < 1/16. Equation (|38|) can be 
rewritten as 



5£ 2/3 = 2/i 2 - h + 2A + 2e 2 + Ah (e - 1) 



p - 4/i G - 4e + 1 



i/p 



p + Ah + At-l 

Using equation (|39|) in the above and retaining terms upto first order in e, we obtain 



(43) 



££ 2/3 = [4h (e - 1) 



p - 4/i D - 4e + 1 
p + 4/i + 4e - 1 



i/p 



(44) 



The two possible values of h a , given by equation (HO), are /i G = 0.25 (1 ± |p|). For 
h = 0.25 (1 — |p|), we choose p > in equation (f44f). This gives 



-|p|e 

-bl 



2|p| -4e 



i/H 



i/H 



(45) 
(46) 



Equation ( IB ) shows that one can satisfy the the condition e — > as £ ^ oo if 1 < |p| 1 
i.e. if 



Vl - 16A < 1 



(47) 



or, in other words, A > 0. (An equivalent result can be obtained for the solution 
h Q = 0.25 (1 + \p\) by choosing p < in equation ([P])). 



Since A > 0, equation ([Tl|) then implies that 



h(h - -) < (48) 

i.e. < /i < 1/2. Thus, the asymptotic values of /i, allowed by the DP solution and 
consistent with the stable clustering hypothesis, lie in the range < h < 1/2. In the 
standard stable clustering scenario, h — > 1 for £ ^> 1. The above analysis clearly rules 
out this value of h as £ —>■ oo. Next, the constraint A > gives 



9QM' 



(49) 



Using equation (|28|) , this gives 



1 + (4 - 7 ) n o - 2S > (50) 

i.e. 



So < ^^n Q + ~ (si) 

We have already shown that S Q > [(4 — 7)/2] II + 1/2 — (3 — 7)/96. The above equation 
(pl|) indicates that only a very narrow range of values is permitted for E Q (in terms of 
n o ), in the non-linear regime. Thus, the DP solution imposes strong constraints on the 
relative values of radial and tangential dispersions; it is not obvious that the dynamics 
of the system will actually cause these constraints to be satisfied. 
Finally, the preceding results do not appear to be influenced by the ansatz h = h(^) 
(although this ansatz is well-motivated, since it appears to be validated in N-body 
simulations; see, for example, [Hamilton et al. 1991| ) as one could have instead carried 
out the analysis using the scaled variable s = xt~ a (with 7 = 2/ (a + 2/3)). In the 
non-linear regime, self-similar evolution implies that h = h(s) and £ = £(s); one can 
clearly write h = h(^) in this regime. 
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4. Conclusions 

In the present work, the BBGKY hierarchy of equations has been revisited and an 
equation derived for the evolution of the dimensionless function h = —(v/ax). The 
assumptions used (following DP) are a hierarchical model for the 3-point correlation 
function and the ansatz that h is a function of the mean 2-point correlation function, 
£, alone, i.e. h = h(C,). No assumption is made regarding the vanishing of velocity 
skewness. In fact, the analysis only uses the zeroth and first moments of the second 
BBGKY equation; it is thus applicable to any form of the skewness (or the higher 
velocity moments) which yields the DP solution for £ and the parallel and perpendicular 
velocity dispersions. The DP similarity solution is then substituted in this equation for 
the h function, in the non-linear regime, and the generalised stable clustering hypothesis 
(h — > const.) used to obtain an expression for the asymptotic value of h, in terms of 7, 
the power law index of clustering and the tangential and radial velocity dispersions. The 
DP solution is found to require that tangential dispersions are larger than radial ones, 
in the strong clustering regime; this can be understood on physical grounds. Finally, 
stability analysis of the solution demonstrates that the allowed asymptotic values of h, 
consistent with the stable clustering hypothesis, lie in the range < h < 1/2. Thus, 
if the DP scale-invariant solution (and the hierarchical model for the 3-pt function) is 
correct, the standard stable clustering picture (h — > 1 as £ — > 00) is not allowed in the 
non-linear regime of structure formation. 
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